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1  INTRODUCTION 

The  following  provides  a  brief  and  informal  account  of  an  analysis  technique  used  by  the  Electronic  Warfare 
Division  { EW' D)  to  produce  high  range  resolution  radar  cross-section  (RCS)  profiles  of  maritime  targets.  VVe 
begin  by  giving  a  simple  theoretical  discussion  of  the  basic  principle  behind  the  technique  employed,  and 
then  move  on  to  describe  the  computer  program  used  to  carry  out  the  first  stage  of  the  RCS  analysis  and 
make  comparison  with  other  worK- 

High  range  resolution  RCS  profiles,  i.e.  RCS  as  a  function  of  range  r,  were  derived  for  a  maritime  target  with 
the  use  of  a  coherent  radar  system  developed  by  the  Microwave  Radar  Division  (MRD).  Synthetic  range 
profiles  were  obtained  for  the  ship  by  exploiting  the  reciprocal  relationship  that  exists,  in  the  context  of  the 
present  work,  between  range  and  frequency'^see  Equation  (3)),  and  by  employing  a  fast  Fourier  transform 
(FFT)  algorithm  to  evaluate  the  inverse  discrete  Fourier  transform  (IDFT).  The  technique  is  referred  to  in 
the  literature  as  stepped  frequency  imag/ng^and  involves  the  recording  of  radar  returns  from  n  pulses,  such 
that  the  carrier  frequency  i/,  for  each  successive  pulse  is  incremented  by  a  fixed  amount  i\i/.  The  tech''tq\ie 
has  been  discussed  by  several  authors,  including  Wehner  et  al.  (1979),  Prickett  and  Chen  (1980).  and  Brvans 
(1986), 


.  y  ,  . 

2  sotAe  theory 

The  purpose  of  this  section  is  to  “fill  in  some  of  the  gaps”  left  by  Wehner  et  al.  (1979).  and  by  Bryans  ( 19861 
in  their  theoretical  discussion  of  stepped  frequency  imaging.  The  treatment  presented  by  the  latter  author 
was  actuallv  based  on  other  work,  including  that  of  Prickett  and  Chen  (1980).  According  to  Bryans  ( 1986  i 
provided  the  duration  of  the  transmitted  microwave  pulse  is  large  enough  to  enclose  the  entire  target,  the 
radar  return  0{k)  say,  from  an  extended  object  at  a  range  R,  can  be  expressed  as  (refer  to  Figure  I  I  .• 

G{k)  =  A  j  F{r]cxp 

where  f  (r)  is  the  targe*  reflectivity  function,  >1  is  a  scale  factor  (a  function  of  power  and  signal  propagation 
for  example),  c  denotes  the  speed  of  light,  and  x  =  1-  The  operating  or  carrier  frequency  i-^.  of  the  radar 

is  given  by 

=  i/Q  -i-  kAt'  (it  =  0,  1, ....  n  -  1).  (2) 

Note  that  i^(n)  is  a  continuous  real  function  of  range  r,  and  G(Jc)  is  a  discrete  complex  function  of  the 
integer  k  (which  in  turn  specifies  the  radar  operating  frequency).  Fuilhetmote,  the  limits  of  the  integral  not 
explicitly  written  in  (I),  are  determined  by  the  width  of  the  transmitted  radar  pulse. 

Wtf  can  also  write  Equation  (1)  as 


+  i:A^)(7?  -  1 


dr, 


F{r)  txp{2'Kxi/Ai)dr, 


(3) 


where  At  =  2(r  -  R)/c,  and  we  have  dropped  the  subscript  notation  on  the  frequency  i/k  That  is.  the 
argument  of  the  exponential  term  in  (1)  and  (3)  represents  the  phase  of  a  returned  signal  relative  to 
^0  =  0.  defined  at  range  r  =  R.  Rearranging  (1),  we  find  that 


where 


//{*)  =  j  F(r)exp  ^  kAv 


)r 


dr, 


4  ir  1 

H(k)  =  A"‘G(il!)exp  1  - — (i/o  -f  kAu)R 
'  c 


(4) 

(5) 
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3  THE  COMPUTER  PROGRAM 


At  the  time  of  writing,  the  most  recent  version  of  the  computer  program  written  to  enable  a  prelinnnary 
analysis  of  the  RCS  data  is  called  HIGH  RES  ANALYSIS  3  (refer  to  the  program  included  as  Appendix  li. 
The  input  data  files  contain  measurements  of  the  in-phasc  (/)  and  quadrature  (Q)  components  of  n  -  26^ 
transmitted  and  received  radar  pulses  per  frame.  The  returned  radar  si^^nrl  C{k),  can  therefore  be  written 
as  0(k)  =  Ik  +  ^Qk-  Header  information  for  each  frame  can  also  be  extracted  from  the  input  files.  An  initial 
radar  operating  frequency  i/©  of  9.16  GHz  was  chosen,  and  a  frequency  increment  of  Ai/  =  1  MHz  used  for 
each  of  the  255  subsequent  pulses.  From  Equation  (10),  the  range  resolution  Ar,  is  then  sci  at  :• 


Ar  = 


3  V  10® 

2  >  256  X  10® 


:::  0-6  m, 


and  yields  an  unambiguous  range  window  (n  -  I)Ar,  of  approximately  150  m,  w-hich  is  large  enough  to 
enclose  the  selected  target. 


The  EFT  routine  used  in  the  program  to  evaluate  Equation  (8)  via 

F(  =  ^  //jt  exp  ^  (/  =  0,  I, .  . ,  n  -  1 ).  <13' 

was  taken  from  p.  754  of  Num^riaii  Recipes  (refer  to  Press  el  al.  1986).  The  storage  convention  ‘.T  the 
input/output  data  array  DATAl.  is  outlined  in  Appendix  11.  Once  the  data  are  windowed  with  the  Hammim: 
window  function,  and  the  n*point  Fourier  transform  applied,  the  amplitude  .4.  associated  with  the  returned 
radar  signal  (a  measure  of  the  power  returned  from  the  target),  is  stored  in  an  output  file,  together  with  the 
appropriate  range  value.  The  above  amplitude  of  the  radar  return  is  presently  defined  as 

.4  =  ;F;i  -  ^^  'Ff  Fi-  (Ml 


That  is,  the  square  root  of  the  sum  of  the  squares  of  the  real  and  imaginary  components  of  the  IDFT  This 
can  later  be  scaled  In  accordance  with  a  standard  calibration  procedure  to  produce  an  estimate  of  kCS  in 
the  units  ofdBm^  for  example. 


The  result  obtained  from  running  the  program  on  an  input  file  is  presented  in  Figure  2.  The  range  profile  wa.« 
plotted  using  the  graphics  package  CRICKET  GRAPH.  Good  agreement  is  found  with  the  corresponding 
plot  for  the  same  data  provided  by  MRD  (see  Figure  3),  although  the  latter  result  was  obtained  with  an 
independent  computer  program  developed  by  that  Division. 


Note  that  the  output  data  files  currently  contain  range  and  radar  return  values  separated  by  a  TAB,  as 
required  by  CRICKET  GRAPH.  The  computer  program  used  in  the  next  stage  of  the  analysis,  i.e.  GRAF^H 
MATCH,  has  been  modified  to  allow  for  this.  Refer  to  Bawden  and  Moon  (1989)  for  further  details  concerning 
the  use  of  GRAPH  MATCH. 


4  CONCLUSIONS 

In  this  research  note  we  have  presented  a  simple  theoretical  discussion  of  a  radar  technique  known  as  stepped 
frequency  imaging.  We  have  demonstrated  the  results  that  one  can  obtain  with  this  technique,  and  shown 
these  results  to  be  in  agreement  with  those  produced  by  MRD  for  a  patlicula'  data  set. 
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Figure  2  Example  of  range  profile  obtained  with  HIGH  RES  ANALYSIS  3  (see  text) 
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Figure  3  The  corresponding  range  profile  obtained  by  MRD  for  Figure  2. 
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APPENDIX  I 

high  rfs  analysis  3 


program  hlgh_re3_analysl3  ; 

{  . 


'  autr.or 

*  dare 

•  pragrara 


P.W.  Taliana  and  s.D.  Elton 
13-MOV-1939 


calculates  the  PCS  proille  ot  a  target  using  a  high 
resolution  radar  technique.  Che  technique  is  based  on 
merhod  and  en-ploys  a  TFT  routine  In  the  analysis 


ISAR 


[  autolinic  runti^ue  units 

(Sr-) 


uses 

Memtypes  ; 

{ ,  qalOcdraw,  oslntC,  tcolint;, 

PasInOut ,  PasConsole,  SAKE,  PasPrinter.  Pac.'ilntC ,  r 

const 

header  start  ■.  "TIME'  ; 
n  -  25o  ; 

?:  -  3.141553  ;■ 

nn3  -  512  ;  {  used  In  ) 

m  *  255  ;  (  j 

Isign  *•  -1  ;  (  routi.ne  ) 


> 


typa 

gldsrray  -  array  tl..nn21  ot  real  ; 
cut_;i:a  -  arraytl. .2,1. .256)  of  real  ; 
ln_file__dati_tytie  »  Record 

header  :  Paciced  Array  (1.  .2561  of  char  ; 
data  ;  array  [1..2561  of 
pac)ced  record 
■r5C_inphase, 

TX~quadrat'are , 

RX'lnphase, 

HXjquadrat'ure  :  byte  ; 
end  ; 

e.hd 


var 

FP.ri.DatAnp 

datal 

MAXJPower 

header 

Ic.l 

lnput_file 

fila_data 

word" 

.number 
in  header 


Out_flle  ; 
gldarray  ; 

Array  ti..2|  of  lo.nginc  ; 
str255  ; 
l.nteger  ; 

file  of  in_flle_data  type  ; 

in  file  dat3_typo  ;  ~ 

STRING;” 

real  ; 

boolean  ; 


Assign  Data 


procedure  assign_data  ; 


I  assign  recieved  signal  ccmoonents 
transformed  ) 


to  a  data  array  that  Is  to  be  Fourier 


begin  (  a33lgn_data  ) 

wrlteln('a33lgn_data')  ; 
for  i  1  to  n  do 
begin 

datal[2-i-i)  :•  fr[2,ll  ;  datil[2*lj  ;•  -fif2,i) 

*  d-  '^'^^'^®^"‘’^'^'>'^*’“>^<3).aat3l(2»l-ll,char01,datalt2*ll)  ;  ) 
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<!r.d  ;  (  J3si.gn_i±5!;a  I 

( . i...  |!a.MMING  . . . . - . ) 

Procedure  H.irrr.ir.g  ; 

{  wiedcw  data  using  Che  Ham.'ing  window  luncrlon  ) 

Var 

a,  1,  Ic  ;  Integer  ; 
wvr.dowl  ;  Real  ; 

Segin  (  Marrrdng  1 
Wciteln  ( '  Ka-nming' )  ; 
a  n  ; 

Tor  1  :•  1  ta  n  do 
Begin 

w-.ndowi  O.r'i-O.-lI  •  cos (2.0  •  pi  •  (i-11  /a)  ; 
darai;;*!-!!  daral  [2*i-ll  •window:.  ; 
daoai;i*i|  datal !2*11 “wlndowl  ; 

Bnd  ; 

ENT)  ;  (  Hirrdng  } 


( . . . .  FAST  FOURIER  TSAKSFCBM  . . •••♦; 

Prooedure  IIt:(vdr  datal:  gldarray;  nn.lsign:  integer)  ; 

{  replaces  IAEA  by  its  discrete  Fourier  trar.slom,  if  ISIGIJ  is  input  as  1;  or 
replaces  IAEA  by  NN  tines  'ts  inverse  dicrete  Fourier  transfom,  i;  I'IGN 
is' .input  as  -1.  DATA  is  a  complex  array  ot  length  NN  or,  equivalently 
a  real  array  oC  length  2*UM.  MN  ausc  be  an  integer  power  o£  2. 

—  Nunerical  Recipes  p  754  ) 
var 

li, Jj,n,nn2.nnax,a, l,lstep,i:  Integer  ; 

wtenp,wr,wpr,wpl.wl,thetai:  double  ;  (  double  precision  for  dhe  trignonetrio  ) 

tenor,  le-Tipi :  real  ;  (  recurrences  1 


becin 

writelni'FFE')  ; 

(  for  1  ;•  1  to  2’nn  do 

begin 

writeln{ 'FFt ' ,  d£tal[il)  ; 
end  :  ) 

nn2  2*nn  : 

1  ;»  1 

for  li  1  to  nn  do 
begin 

1  2*11-1  ; 

if  (j>l)  then  begin 
tempr  datal [ 31  ; 
tempi  :«  datal [j*!!  : 
datal [ 31  :■  datal (11  ; 
datal(3+ll  datal(l+ll  ; 

datal (il  :•  tempt: 
datal[l+ll  tempi 

end  ; 

m  nn2  dlv  2  ; 

while  i  (ra  >-  2)  and  ( j  >  m) )  do 
begin 

3  3-in  ; 
ra  ;•  ra  div  2 
end  ; 

3  3*ra 
end  ; 

mmax  :■  2  ; 

while  (nn2  >  ramax)  do 


(  t.his  is  the  bit-reversal  section  of  1 
(  the  routine  1 


(  exchange  the  complex  numbers  1 


(  here  begins  the  Dan.elson-Lanctos  section  of  the  r 
(  outer  looo  executed  log(2)  NN  times  1 
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ure  Write  lata  I'at  ; 


freqir.c  = 


ir.sax  :»  V  ; 

;  :■  -r,  dlv  2  to  n  iiv  2  -  1  do 
becin 

range  :»  (C'lndex)  /  (2.0’freqir.c*n)  ;  {  lonver 

index  index  *  1  ;  (  “FT 

if  :  <  0  tnen  i  :»  :  •  n 
else  i  ;•  3  ; 

a-Tp  :■  dAtai[2*'i-i;'*datai!2*l-ll-datal.'2*iI*data 
ar.pl  :■  10.0  •  In  lartp)  •  0.A043  ; 
anp2  r-  round(ar.pll 

{  writeln  (oucput_file,  range:  4:  l^c.tar  (9)  .arp2: 3)  ;  ) 
writein  (cutput^^flle/  range :  5:2, o^.ar  (9)  ,a.ttpl :  5 : 2)  • 
end  ;  " 

writeln  (output  file,  {2C*ln(DatJ«TOi; l.i; )  •0.4343) :  6,crir  (?) . 

(20'ln(DatArro'i2,  i; )  •0.4342)  :6,Chr  (9) 

)  ;) 


d  index  associated  witn 
to  dC'ilfied  range  • 


i:2*i:  ; 

i  log(xj  •  ln(x» • : . 434  2) 


close )output_flle)  ; 
e.nd  f  write  data  out  ) 


Main 


Progam 


Pegin 

reset (tnput^flie, 'rcsl00.dat' )  : 
for  k  :■  i  to  2  do 
term 

Read_2ata_:n(  k  )  ; 
assign_data  ; 

Haiming  ; 

“FT (datal, nn, Isign)  ; 
Write__lata_out  ; 
end  ; 
readln 

end. 
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APPENDIX  II 

FFT  STORAGE  CONVENTION 


Figure  II  I  illustraies  ihe  storage  convenlion  for  ihe  inpiil,  output  data  array  DATAl  contained  in  the 
program  HIGH  RES  ANALYSIS  3.  The  input  array  {refer  to  (a|)  contains  N  (a  power  of  2)  complex  data 
samples  in  a  real  array  of  length  2.V.  will,  real  and  imaginary  components  alternating.  The  output  arrav 
(refer  to  (b))  contains  the  complex  Fourier  spectrum.  Again  the  real  and  imaginary  components  alternate. 
The  array  begins  with  zero  frequency,  and  builds  up  to  the  most  positive  frequency.  Negative  frequencies 
then  follow,  starting  with  the  second  most  negative  frequency  (remember  the  most  positive  and  negative 
frequencies  are  ambiguous),  up  to  the  frequency  just  below  zero. 
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Figure  II.l  Storage  convention  for  input/output  data  array  DATAl  (after  Press  et  al.  1986). 


ERL-0502-RN 


DISTRIBUTION 

DEPARTMENT  OF  DEFENCE 

Defence  Science  And  Technology  OrgAnisation 
Chief  Defence  Scientist  1 

First  Assistant  Secretary  Science  Policy  J 
Counsellor  Defence  Science,  London 
Counsellor  Defence  Science,  Washington  j 
Electronics  Research  Laboratory 

Director,  Electronics  Research  Laboratory 

Chief,  Electronic  Warfare  Division 

Research  Leader,  Electronic  Countermeasures 

Head,  Electronic  Countermeasures  Group 

Dr  P.R.  Gerhardy,  Electronic  Countermeasures  Group 

Dr  T.T.  Moon,  Electronic  Countermeasures  Group 

Mr  P.W.  Taliana,  Electronic  Countermeasures  Group 

Surveillance  Research  Laboratory 
Chief,  Microwave  Radar  Division 
Research  Leader,  Microwave  Radar 
Head,  Electromagnetics  Group 
Dr  W.C.  Anderson,  Electromagnetics  Group 
Dt  L.  Berzins,  Radar  Techniques  Group 

Weapons  Systems  Research  Laboratory 
Director,  Weapons  Systems  Research  Laboratory 
Chief,  Guided  Weapons  Division 

Joint  Intelligence  Organisation  (DSTI) 

Air  Office 

Air  Force  Scientific  Adviser 
Navy  Office 

Navy  Scientific  Adviser 
Army  Office 

Army  Scientific  Adviser 

Libraries  and  Information  Services 
Defence  Library,  Campbell  Park 

Document  Exchange  Centre,  Defence  Information  Services  and 
Science  Liaison  Branch  for  Microfilming 
National  Library  of  Australia 
United  States  Defence  Technical  Information  Centre 
United  Kingdom  Defence  Reseaich  Information  Centre 
Director,  Scientific  Information  Services,  Canada 
New  Zealand  Ministry  of  Defence 
British  Library,  Document  Supply  Centre  (UK) 

Institution  of  Electrical  Engineers  (UK) 

Australian  Defence  Force  Academy  Library 

Main  Library,  Defence  Science  4r  Technology  Organisation  Salisbury 
Library,  Aeronautical  Research  Laboratories 


Copy  No. 


1 

Control 
Sheet  Only 

2 

3 

4 

5 

8 


9 

10 

11 

12 

13 


14 

15 

16 

Control 
Sheet  Only 

Control 
Sheet  Only 

Control 
Sheet  Only 

17 

18 
19 

20-31 

32-A^ 

44 

45 

46 

47 
46 

49-50 

51 


ERL-0502-RN 


Library,  Materials  Research  Laboratories 
Librarian,  Defence  Signals  Directorate 


DOCUMENT  CONTROL  DATASHEET 


Security  classification  of  this  page  ;  UNCLASSIFIED 


I  DOCUMENT  NUMBEFIS 
AR 

Number :  AR-0()5-982 

Series 

Number :  ERL-0502-RN 


2  SECURITY  CUSSIFICATION 

a.  Complete 

Document  :  Unclassified 

b.  Title  in 

Isolation ;  Unclassified 


c.  Summary  in 

Isolation :  Unclassified 


Other 
Numbers  : 


2  DOWNGRADING /  DELIMITING  INSTRUCTIONS 


TITLE  .  ’  '  . 

HIGH  RANGE  RESOLUTION  RADAR  CROSS-SECTION  PROFILES  -  A  PRELIMINARY 
ANALYSIS  (U) 


5  PERSONAL  AUTHOR  (Si 


®  '  DOCUMENT  DATE 
January  1990 


Stephen  D.  Elton 


I 


2  I  7.  1  TOTAL  NUMBER 
OF  PAGES 


7.  2  NUMBER  OF 

REFERENCES  5 


^  8.  1  CORPORATE  AUTHOR  (S) 
Electronics  Research  Laboratory 


^  n^FERENCE  NUMBERS 
I  a.  Task : 


b.  Sponsoring  Agency  : 


8.  2  DOCUMENT  SERIES 
and  NUMBER 
Research  Note 
0502 


;  COST  CODE 

I 

i  _ _ 


IMPRINT  (Publishing  organisation)  i 

— - - 1 

Defence  Science  and  Technology  I 

Organisation  Salisbury  j 


I  COMPUTER  PROGRAM  (S) 
(Title  (s)  and  language  (s)) 


!  RELEASE  LIMITATIONS  (of  the  document) 
Approved  for  Public  Release. 


Security  classification  of  this  page  : 


UNCLASSIFIED 


Security  classification  of  this  page  :  UNCLASSIFIED 

ANNOUNCEMENT  LIMITATIONS  (of  the  information  on  these  pages) 


No  limitation 


^  DESCRIPTORS  ' 
a.  EJC  Thesaurus 

Terms  Radar  cross  sections 


COSATI  CODES 

1709 


b.  Non  ■  Thesaurus 
Terms 


Stepped  frequency  imaging 


SUMMARY  OR  ABSTRACT 

(if  this  is  security  classified,  the  announcement  of  this  report  will  be  similarly  classified) 


We  briefly  review  the  technique  of  stepped  frequenq/  imaging,  and  describe  the  implementation  of  a 
computer  program  that  calculates  high  range  resolution  radar  cross-section  profiles  of  maritime 
targets. 


Security  classification  of  this  page  : 


UNCLASSIFIED 


